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Abstract. 

Measurements of the Sunyaev-Zeldovich (hereafter SZ) dis- 
torsion of the cosmic microwave background can give inter- 
esting physical informations on clusters of galaxies, provided 
that the electronic temperature of the gas is known. Previous 
attempts to do so have used the electronic temperature de- 
termination obtained from the X-ray spectra. However, if the 
intergalactic gas is not homogeneous, the X-ray emission will 
trace the denser component, and the temperature determina- 
tion may not be relevant for the lower density gas which is 
dominating the SZ measurements. Moreover, the X-ray bright- 
ness decreases very rapidly with the distance, which is not the 
case for the SZ effect. Distant clusters might be detected from 
SZ measurements, whereas they are inaccessible to X-ray ob- 
servations. For these reasons, we have investigated the possibil- 
ity to derive the electronic temperature of the gas from the SZ 
measurements in the submiUimeter range (A ~ 300 — 600^m) . 
We show that given the sensitivities of the future submiUime- 
ter space missions Planck Surveyor and FIRST, the electronic 
temperature of massive clusters {Ycenter = 3 x 10"*) can be de- 
termined with an accuracy ranging from 1 to 4 keV depending 
on its distance and the data available. 
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1. Introduction 

The intergalactic medium is a stro ng source of diffuse X - 
ray radiation by free-free emission (Jones & Forman 1984) 



It is observed at submiUimeter wavelengths too, via the 
Sunyaev-Zeldovich (hereafter SZ) effect: a spectral distorsion 
of the Cosmic Microwave Background (hereafter CMB) due 
to the interaction of the electr ons of the hot ionized gas 
with the photons of the CMB dZel'dovich fc Sunyaev 1969| , 



Sunyaev & Zeldovich 1972). If the electronic temperature is 
determined from X-ray data, together with a model of the 
gas distribution, SZ data allow to derive the Compton opti- 
cal depth (r) of the intergalactic gas. This parameter directly 
provides the gas mass, if it is integrated over solid angles. The 
association of X-ray and SZ data allows to estimate the Hub- 
ble constant. Ho, independently of the usual standard candles 



methods (see Holzapfel et al. 1997 for instance). The peculiar 
velocity of several clusters can also be derived from the Doppler 
effect, so that it should be possible to detect the large scale 
gravitational field which is produced by the dark matter. 

An important feature of the SZ brightness is that it is an 
absorption effect on the CMB, which intensity is independent 
of z (if no evolution of the cluster is assumed). On the contrary, 
the cluster's X-ray surface brightness decreases with respect 
to the usual (1 -|- z)~^ expansion factor and with respect to 
an additionnal exponential factor, exp{—E{l + z)/kTE), which 
becomes important for very distant clusters. 

In this paper, we investigate how the gas cluster tempera- 
ture can be recovered from the SZ measurements themselves. 
In section 2, we present a simple Monte-Carlo method which 
allows to obtain the exact shape of SZ spectra, taking into 
account the temperature dependence which is ignored in the 
usual analytical expression. We emphasize the spectral effect 
due to the gas temperature which shows up in the submiUime- 
ter domain. In section 3 we quantify the error on the deter- 
mination of the gas temperature using SZ measurements only, 
considering limitations of the instruments sensitivities and of 
the various foreground and background emissions. 

2. Exact calculations of the SZ effect 



The frequency dependency of the SZ effect is usually ap- 
proximated by a solution of the Kompaneets equation 
(Kompaneets 1972) whic h is a second order approx imation of 
the Boltzmann equation ( Rybicki fc Lightman 1979 ). This so- 
lution needs: 
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for the kinetic effect, 
f{x) and a{x) are analytic functions of the dimensionless fre- 
quency: 
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y is the Comptonization parameter of the cluster: 



y 



kT, 



(4') 3.1. Precision of X-ray determinations 



Te is the electronic temperature of the intergalactic gas, 
Vp is the pecuhar velocity of the cluster, r is the Compton 
optical depth of the intergalactic gas. T,. = 2.726K is the CMB 
temperature. 

Photons of the CMB are scattered from low frequencies to 
higher frequencies. The cross-over frequency is around v = 217 
GHz (e.g. A = 1.38 mm). 

In case of millimeter and submillimeter SZ data, the ana- 
lytic approximation yield to analysis errors. Indeed, the elec- 
tronic temperature level implies weakly relativistic velocities 
for the electrons. A few auth ors have work e d on the rela tivistic 
corrections of the SZ effe ct ([Wright 1979|, [Fabbri IQSi]). Most 



recently, Rephaeli ( 1995 ) has compiled previous works to de- 
velop a semi-analytical treatm ent of the SZ effect a nd Ste bbins 
( [1997| ), ChaUinor & Lassenby ( [l997| ) and Itoh et al. ( [19971 ) have 
extended analytically the Kompaneets equation. However, the 
exact SZ spectra can be obtained using a simple Monte-Carlo 
method which numerically integrates the transfer equation: 



^ = n. fdp j dnpm^[h{u^)~ I..{u))] 



(5) 



where Pe(/3) is the velocity distribution of the electrons, 

13=—, "Tq is the differential scattering cross section, Ii, is the 
intensity of the radiation at the frequency v (before scattering) 
and v\ (after scattering), n^ is the density of the electrons. 
The frequency shift of the photons is given by: 
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where B = (3 cos 6 and Bi = j3 cos 9i and cos a = cos 6 cos 6i + 
sinSsin^i cos{4> — 4>\), {0,4') ^^d {9i,<j)\) are the angles between 
the electron's and the photon's directions of propagation, re- 
spectively before and after the scattering. 

The exact differential cross-se ction has to be used 



( Podzniakov, Sobol fc Sunyaev 1983 ) 
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Te = e^ / rUeC? is the classical radius of the electron and 
7 = (1 + /3^)~^ is the Lorentz's factor. 

As shown on Fig. [l[, the spectra that we obtain ed ar e 
in accordance with those calculated by Rephaeli (1995). 
Those spectra can be obtained from our anonymous ftp site 
(ftp.cesr.fr/pub/astrophysique/sz/). We will show hereafter 
that the SZ dependency on Te can be used to derive the in- 
tracluster gas temperature from submillimeter data (see Fig. 

3. Determinations of the electronic temperature 

In a few years the Planck Surveyor and FIRST missions will 
detect the SZ effect with a very good sensitivity. It is likely that 
new distant clusters will be identified with SZ data and that 
these clusters will remain out of reach for the X-ray and/or 
optical telescopes. In the following, we discuss the two ways of 
determining the electron temperature for galaxy clusters from 
X-ray or SZ data. 



Using the XSPEC software (Arnaud 1996) we simulated a 25 



kilosecond ASCA (GIS2) observation of a galaxy cluster with 
the following spectral parameters: Bremsstrahlung tempera- 
ture (Te) of 8 keV, an unabsorbed 1-10 keV flux of 6 x 10"" 
ergs s~^ cm~^ corresponding to a rich cluster with Q < z < 0.2 
observed through a column density of 10^^ H atoms cm"^ por 
such an observation, Te is recovered with an error of typi- 
cally 0.5 keV at 68% confidence level. These errors are consis- 
tent with those derived from actual observations (David 1993). 
Given the ASCA energy range (0.1-10 keV), the error on Te 
should increase with Te . This is simply because as Te increases 
the cutoff in the Bremsstrahlung spectrum moves at higher en- 
ergies, and approaches or even exceeds the high energy thresh- 
old of the instrument where its sensitivity drops sharply. Con- 
sequently, fitting the spectrum will tend to underestimate Te. 
This is a major limiting factor which is primarily related to 
the energy coverage of the instrument. This will also apply to 
future instruments, like XMM, although their sensitivities are 
much better than ASCA. For instance simulating an ASCA 
observation with the same input spectrum as above but with 
Te = 12 keV, the fit recovers Te as 10.9 ± 0.9 keV (90% con- 
fidence) barely consistent with the input value. Obviously the 
recovered Te and associated errors via X-ray observations de- 
pends also on the input flux which has to be compared with 
the instrument detection sensitivity. As the X-ray flux at earth 
decreases sensitively with the redshift ( e.g. z>l), for distant 
clusters the accuracy in determining Te will also decrease. For 
instance, Hattori et al. ( 1997 ) fitted the ASCA data of the 
AXJ2019+1127 X-ray cluster to z = 0.94 and Te = %.&ttl'^- 
In the case of an inhomogeneous intracluster gas, if we con- 
sidered the temperature-density correlation, a Te determina- 
tion via X-ray emission (oc J n^^ipjdl) may miss a low density 
component of the gas, so that an independent temperature de- 
termination, via SZ measurements (mostly sensitive to a low 
component, e.g. oc f ne{l)dl) is also interesting. 

3.2. Temperature determination with Planck and FIRST 

In the following, we have simulated SZ measurements of a rich 
cluster performed with Planck and FIRST. The photometric 
bands used are those of the HFI for Planck, PHOC and SPIRE 
for FIRST. Their characteristics are summarized in Table m. 
They are taken from the COBRAS/SAMBA Phase A report 
concerning Planck (except for the additional bolometer chan- 
nel at 100 GHz). Concerning FIRST, we have assumed that 
the noise level used is the quadratic sum of the photon noise, 
NEPphot = hv{^ Au€rin{l + erju))^^'^ and the detector noise 
NEPdet = 3 X 10"" W/VHz. Sn/X^ = 1 at the diffraction 
limit, IS.V jv ~ 0.3, e = 0.03 is the telescope emissivity, r] = 0.5 
is the system response (including transmission and detector ef- 

ficiency), n = 2/(e '''^t^i — 1) is the photon phase space occu- 
pation number for the telescope thermal emission, Tt^i = 80 K 
(sateUite at I/2 Sun-Earth Lagrangian point.). The wavelengths 
of the FIRST lower frequency channels have been adjusted to 
the Planck higher frequency bands. 

We estimate the precision of the temperature determina- 
tion by repetitive least square fits on simulated submillimeter 
and millimeter data: Planck alone or Planck plus FIRST com- 
bination. We take into account both the instrumental noise and 
the contaminating sky emissions: galactic dust, free-free, syn- 
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Table 1. Characteristics of the Planck Surveyor's and the FIRST'S passbands used. The Noise Equivalent Brightness (e.g. NEB) 
is given for the nominal beams for Planck and for a virtual beam for FIRST (i'. The rms contribution of the backgrounds and 
foregrounds are given for the beam used (see caption of Fig. 1 and text). Dilution factors applied to Ycenter = 3 x 10"* are given 
for redshifts 2: = 0.1 and z = 1 and for a cluster with Re — 0.3 Mpc and fitot ~ 0.3, Ho = 70 km/Mpc/s. 

( ' A virtual beam of 4.5 arcmin (e.g. the Planck submillimeter beam) is used in simulations. It will be obtained by summation 
of adjacent sky pixels. The NEB in each bands is the sum of photon noise and detector noise, as determined in section 3.2.1, 
modified by a factor equal to the ratio of the diffraction limit beam to virtual beam diameters 



chrotron, Ultra-Luminous InfraRed Galaxies (ULIRGs), and 
CMB. We did not attempt to determine the velocity through 
the kinetic effect since it is spectraly identical to the pri- 
mary CMB and has alre ady been extensively stu died (see 



Haehnelt fc Tegmark 1996 and Aghanim et al. 1997) 



3.2.1. Simulation of clusters observations 

We simulate a massive cluster, Ycenter = 3 x 10"*, with isother- 
mal /3 density gas profile, Te = 8keV, /3 = 2/3, Re = 0.3 Mpc, 
no = 1.9 X 10~^cm~"^. We take into account the dilution of the 
cluster SZ emission in the instrumental beam by proper inte- 
gration of the SZ profile over the beam (we assume the density 
cutoff at 15 Re)- Concerning the FIRST channels we use with 
a virtual beam equal to the Planck submillimeter beam (4.5 
arcmin) which will be obtained by the summation of adjacent 
sky pixels. The surface brightness noise level (NEB) is cor- 
respondingly improved from the diffraction limit value by a 
factor equal to the ratio of the diffraction limit beam to the 
virtual beam. The signal is then correctly integrated over the 
finite bandpass of the instruments. 

We then add the following astrophysical contributions: 

- A dust signal with rms level lioon =0.3 MJy/sr, value 
which is not exceeded on 19% of the whole sky(This percentage 
is determined from an all sky map of the rms 100 fim fluctu- 
ations calculated in bins of 0.7° x 0.7° from the ISSA IRAS 
maps). A spectral characteristics T — 1 7.5K , n — 2, typical 
of high latitude clouds (Boulanger et al. 199q). The variability 
of the dust spectrum is taken into account by introducing a 
random fluctuation in the spectral index. This is fixed to 10% 
rms of the average value (n = 2), in accordance with mea- 
surements of the 2 meter submillimeter stratospheric telescope 
PRONAOS on high latitude galactic cirrus (Bernard et al., in 
preparation). 

- Fluctuations of the integrated contribution from the back- 
ground ULIRGs. Their rms l evel is i ntegrated from the number 



counts of Guiderdoni et al. 



(1997) 



The fiux limit is set such 



that the probability to find a source brighter than that limit 
within the beam is smaller than 10% for the Planck's bands. 
For FIRST, considering the high angular resolution available, 
we assume that all sources above the 3a noise level are detected 



and subtracted before summation in the 4.5 arcmin virtual 
beam. 

- Fluctuations of the primary CMB at a rms level of 
AT/T = 30 x 10"*^ (e.g. A^p ~ 500 km/s for our cluster). 

- free-free and synchrotron: Their level and spectral be- 
haviour are taken from the COBE determination (Bennett et 
al. 1992): 7/iii' rms free- free at 53 GHz with spectral index n 



= -0.16, 4fiK rms synchrotron at 53 GHz with spectral index 
n — -0.9. (The COBE measurement h as been extrapolated t o 
sub-degree angular scales assuming l~^, Bersanelli et al. 1996| )^ 

A random instrumental noise with rms level as in Table hi 
is then added to the astrophysical signal. Finally, a 1% relative 
error is randomly added to all bands to take into account the 
band to band calibration uncertainty. 

The different emissions and noise levels are shown in Fig. 
H together with the instrumental band positions. 

3.2.2. Restoration of the cluster's parameters 

We fit a five parameters model to the simulated multi-band 
measurements. This model is the sum of : SZ thermal effect 
(Te, Fcenter), primary CMB (AT/T) and dust {Cd,n). Ca and 
n are respectively the level and spectral index of the dust. 
The model also take into account the integration over different 
spectral bands. 

The errors on the restored parameters are estimated from 
the statistics of the multiple simulations and fits. The ULIRGs 
are not explicitly included in the separation because their spec- 
trum is very close to galactic dust spectrum, however they are 
partly taken into account through the variable spectral index 
of the fitted dust. 

4. Results and discussion 

Due to the likely curvature of the universe, a massive cluster 
as the one considered in the previous section {r^ > 0.3 Mpc), 
will never show a SZ profile having a width at half maximum 
less than 2 arcmin. The integration time (~ 10 s) is the main 
limiting factor concerning the use of the Planck survey to mea- 
sure the SZ effect with a high accuracy. Actually, this should 



be overcome by the FIRST telescope. Whatever are the mis- 
sions final concepts and observing strategies, data from the 
two satellites will be available in a few years, and will provide 
a unique tool to study distant clusters of galaxies via precise 
measurements of the submillimeter SZ effect. 

The accuracy on the determination of Te is strongly depen- 
dent on the integra tion time. The FIRST Deep survey, a fore- 
seen key program (Beckwith et al. 1993), is expected to cover 
100 square degrees on the sky, for a total time of 2 x 10^ seconds 
that corresponds to an integration time of 500s per position, 
which is the value that we will use in the following. According 
to the Planck's mission observing strategy, some parts of the 
sky will be observed with an integration time a factor of 10 
higher than the nominal time, (about 100 squared degrees on 
the sky.). 

We restore the gas temperature for these two integration 
times in the cases of Planck alone and the FIRST/Planck data 
combination. For Planck alone, effective beams are all fixed to 
10.6 arcmin (the 3mm beam of Planck) to obtained a better 
sensitivity. 

Results are summarized in Table H. 

The accuracy on the estimation of the SZ parameters and 
the number of clusters that will be detectable depends on the 
observing strategy. A first one would be to search for rich clus- 
ters in the FIRST deep survey fields, completed by the corre- 
sponding positions of the Planck survey. One can expect that 
part of the FIRST deep survey will be performed in a region 
where Planck has a higher integration time per pixel (P2/F 
case in Table 2) . The precision on the temperature determina- 
tion is in this case of 1.2 keV {z — 0.1) to 3.9 {z = 1.). Massive 
distant clusters can be identified in the FIRST data as positive 
extended excesses at the larger wavelengths, after removal of 
the dust component. (The dust contamination should not be a 
problem since the deep survey field will be selected among the 
best galactic windows). 

A second strategy will consist of a search for rich clusters 
in the Planck survey, completed by pointed observations with 
FIRST. The distant rich clusters will show almost no spatial 
extension in the Planck survey {FWHAlsz — S.56c — 2.2 ar- 
cmin at 2 = 1) and will be recognized among the faint clusters 
detected in the Planck survey, because they will have no or 
only faint X-ray counterparts. Actually, massive clusters (e.g: 
2/ > 3. X 10~*) will be detectable in the Planck survey data 
at all redshifts; the most extreme dilution factor for such a 
cluster is about 0.18 (10.6 arcmin beam) and the instrument 
sensitivity to SZ effect is of order 1.34 x 10~® (no dilution). For 
each cluster, depending on the detector array size, a five ar- 
cmin field on the cluster (plus comparison field) can be covered 
with FIRST within a few hours of observation. The precision 
on the temperature determination will be of order 1.2 to 1.4 
keV for a cluster ai z — 0.1, depending on the Planck integra- 
tion time on this cluster. For very distant clusters the dilution 
is important and the precision is limited to ~4 keV due to 
the contamination of the infrared background galaxies. This 
result, which is comparable to th e prec ision of the current X- 
ray determinations (Hattori et al. 1997), could be improved by 



cleaning properly the emission in the Planck bandpass from 
the infrared galaxies contamination. 

Depending on the age of the universe, its geometry and the 
epoch of cluster formation, th ere may be a few or a large num- 
ber of distant rich clusters (Oukbir & Blanchard 1997). Al- 
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Table 2. 1 a statistical error (keV) on the determination 
of the gas temperature (8 keV cluster) for different config- 
urations. For Planck alone, PI: nominal integration time. 
P2:10xnominal integration time. For Planck/FIRST, Pl/F: 
Planck's nominal integration time and First deep survey (DS) 
for FIRST. P2/F: Planck 10 x nominal time and First deep sur- 
vey 



though cosmological standard models (CDM and f2 ~ 1) favour 
the formation of rich structures only at a recent time, and seem 
t o be supported by the X-ray cluster's distrib ution at z < 0.5 
(Bartlett & Silk 1994, Luppino & Gioia 1995), several detec- 



tions of galaxies and structures have been reported at very high 
(I < z < 5): 



Le Fevre et al. (1996 



Pascarelle et al. (19961 ^ 
and Pey et al 



Malkan et al. (1996) 



19W 



at optical wave- 
leng ths, t he X-ray detection of AXJ2019-I-1127 by Hattori et 
al. (1997 ) at z=0.94. Most une xpected have been the detec- 



tions by lones et al. (1997) and Richards et al. (1997) of neg- 
ative decrements at centimeter wavelengths in the direction of 
distant quasars. These decrements can be attributed to massive 
clusters at z > 1, u nobse rvable at other wavelengths. 

Korolyov et al. ( 1986 ) were the first to write that it is easier 
to detect high redshift clusters of galaxies in microwave spectral 
band than in X-rays. If this is confirmed, it means that a fairly 
large number of rich clusters may be detectable in the Planck 
survey and selected on the basis of low dust contamination and 
faint (or no) X ray counterparts for observations with FIRST. 
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Fig. 1. Comparaison of the SZ spectra obtained using the ana- 
lytical approximation (dotted lines) and those obtained using 
the exact Monte Carlo calculation (solid lines). Spectra are 
plotted for clusters with 5, 10, 15 and 20 keV temperature (in 
units {hcf/[2{kTef]). 
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Fig. 2. The temperature effect in the thermal SZ effect 
(ASZ^difference between exact and approximation calcula- 
tion): solid line for the positive part and dotted-dashed line for 
the negative part. The exact SZ effect (long dash line) and the 
analytical approximation (dotted line) are overplotted. ( Te = 8 
keV, Ycenter = 3 X 10"^). Other contributions have been plot- 
ted: The CMB at the level of AT/T = 30 x 10"", the galac- 
tic free-free plus synchrotron emission respectively normalized 
at the level of 7fiK and 4^K at 53 GHz (dashed lines). The 
galactic dust (with Td = 17.5 K , n = 2 and normalized to 
0.3 MJy/sr at 100 /im) and the background infrared galaxies 
beam to beam fluctuations (3 dot-dashed lines). (The back- 
ground galaxies fluctuations are plotted as two curves, for the 
Planck and for the FIRST channels, as described in the text.). 
Horizontal bars represent the Plank and the FIRST sensitivi- 
ties levels. (Dashed bars for the high sensitivity regions of the 
Planck survey.) 
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